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ABSTRACT 

Aims. Using a counter-dispersed slitless spectroscopy technique, we detect and measure the line-of-sight velocities of 187 planetary 
nebulae (PNe) around one of the nearest cD galaxies, NGC 1399, with FORS1 on the VLT. 

Methods. We describe the method for identifying and classifying the emission-line sources and the procedure for computing their 
J2000 coordinates and velocities. The number of PN detections and the errors in the velocity measurements (37 km s -1 ) indicate that 
this technique is comparable to other methods, such as that described by Teodorescu et al. (2005). 

Results. We present the spatial distribution of the PNe and a basic analysis of their velocities. The PN two-dimensional veloc- 
ity field shows marginal rotation consistent with other studies. We also find a low-velocity substructure in the halo and a flatter 
velocity-dispersion profile compared to previous observations that extends to ~400". The detection of a low-velocity subcomponent 
underscores the importance of discrete velocity tracers for the detection of un-mixed components. 

Conclusions. The new velocity-dispersion profile is in good agreement with revised velocity dispersions for the red globular clusters 
in NGC 1399, using the data of Schuberth et al. (2009). The outer parts of this profile are consistent with one of the dynamical models 
of Ki'onawitter et al. (2000), which corresponds to a circular velocity of ^340 km s _I and a rescaled B-band mass-to-light ratio of -20 
at T radius. These measurements trace the kinematics of the outer halo and disentangle the heterogenous populations in the Fornax 
Cluster core. The new data set the stage for a revised dynamical model of the outer halo of NGC 1399. 

Key words. Galaxies: individual: NGC1399 - Galaxies: elliptical and lenticular, cD - Techniques: spectroscopic 



1. Introduction 

An important key to understanding the formation and evolu- 
tion of galaxies is to examine the kinematics of the halo where 
the orbits have long dynamical timescales and the mass pro- 
file is dominated by dark matter. In disk galaxies, HI measure- 
ments generally provide line-of-s ight velocities out to the edges 
dRubin et alj|1977t lBosmalll978l) . Measuring kinematics in the 
outer parts of early-type galaxies is more complicated since HI 
measurements are usually unavailable, and the signal-to-noise 
of the integrated light spectra drops off before the mass pro- 
file flattens. Specifi cally, long-slit spectra pro vide measurements 
for the first 2-3 R e (iKronawitter et al.ll2000t) . multi-object spec- 
tros copy provides line -of-sight velocity distributions out to 3 
R e (Proctor et al. I2009h and IFUs can take measurements at 4 
R e (iMurphv & Gebhardjl2007t IWeiimans et al.ll2009l) . Discrete 
tracers like planetary nebula provide the necessary kinematic in- 
formation to examine the regions beyond these limits and learn 
about th e formation his t ory of the outer regi ons of elliptical 
galaxies dHui et al.lll99llArnaboldi et al.ll 1996b . 

The planetary nebula (PN) stage occurs at the end of a 
star's life when the atmosphere has been ejected and is subse- 
quently ionized by UV radiation from the exposed stellar core. 
Due to the extreme rarefaction of the gas in the atmosphere, 
these objects emit up to 15% of their light at [OIII] (/15007) 
dDopita et al.l 1 1992b . which allows their detection and veloc- 



ity m easurement out to a distance of 100 Mpc dGerhard et al.l 
12005b . Moreover, t he PN number dens ity traces the light in 
early-type galaxies dCoccato et al.ll2009l) . Because they follow 
the same spatial distribution as stars, PNe provide kinematics 
that are directly related to integrated light measurements. We 
are able to use planetary nebulae as kinematic tracers out to 
more than 7 R e for measuring the dark m atter distribution and 
angular momentum content in the halo ( Mendez et alj 2001 



* Based on observations made with the VLT at Paranal Observatory 
under programme ID 066. B -027 8 



i Romanowskv et al.ll2.Q03t iDouglas et alj|2007t Ide Lorenzi et al.l 
l2008ll2009LlNapolitano et al.ll2009UCoccato et alJl2009h . 

In general, globular clusters follow a different spatial distri- 
bution so their kinematics do not trace the same population as 
the st ars (lHarrislll99H I Ashman & Zepfll 19981 iBrodie & StraderJ 
12006b . The globular cluster system is dominated by the blue 
component of a bimodal color distribution in most galaxies. The 
radial distribution of this component is flatter than the integrated 
light from the stars. However, in many galaxies including NGC 
1399, the red globular cluster sample does follow th e integrated 
light dBrodie & StradetteOO&rSchuberfh et al.l2009b . and we ex- 
pect the kinematics of this subset to be comparable to those of 
the integrated light and the PNe. 

Traditi onally, PNe are detec ted using an on-band/off-band 
technique dCiardullo et al.|[l989h . The sky is imaged through a 
narrow-band filter centered at the redshifted [OIII] A 5007 line 
and then again in a broad band that does not contain any bright 
emission lines. PNe are too faint to detect in the broad-band im- 
age, but they appear in the narrow -band image. Foreground stars 
appear in both narrow- and broad-band images, so the wide fil- 
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Fig. 1. DSS image of NGC 1399 overlaid with the five observed 
fields. The fields were intentionally overlapped in order to con- 
firm consistent astrometry and velocity calibration across the ob- 
servations. NGC 1404 is seen in field 7. 



ter is used to distinguish them from the PNe. Once the posi- 
tions are measured from the on-band image, the velocities are 
then typically obtained using multi-slit o r fibre spectroscopy 
dArnaboldi et al.ll2004t iDohertv et alj|2009l) . 

More recently, a number of techniques have been developed 
that allow detection an d velocity measurem ent to be done in 
the same observing run. Mendez et al.l d2001l) identify PNe us- 
ing narrow- and broad-band images, but then follow with a dis- 
persed image in the same run. The velocity of each PN can be 
calculated from the difference between the object's position in 
the narrow-band frame and the dispersed frame. This technique 
has also been applied to NGC 1344, another Fornax galaxy, by 
iTeodorescu et al.ld2005f) . We compare the efficiency of those ob- 
servations to ours in 33.3.11 

The PN Spectrograph (PN.S) is an instrument on the William 
Herschel Telesc ope dedicated to m easuring PN velocities in 
nearby galaxies dDouglas et al.ll2002h . The pupil is split in half 
before being dispersed in opposite directions in twin spectro- 
graphs. A combination of the two exposures, which are simulta- 
neously generated, allows identification of emission-line objects 
and the measurement of their velocities using the separation be- 
tween the positions in the two spectral images. 

For extremely faint samples in distant galaxies, observing 
the entire field through slits significantly reduces the signal from 
the sky. Multislit Imaging Spectroscopy involves a grid of slits 
that are stepped until the entire field has been spectrally imaged. 
These types of observatio ns have detected samples of intraclus- 
ter PNe 100 Mpc away dGerhard et al.l 12005b lArnaboldi et all 
l2007tlVentimidia et al.ll2008l) . 

The counter-dispersed slitless technique in this paper uses 
only two exposures for each field to obtain positions and veloc- 
ities of PNe in the cD galaxy, NGC 1399. The field is first ob- 
served with a dispersed image. Next the spectrograph is rotated 
180 degrees and the same field is exposed again, this time with 
the dispersion in the opposition direction. As in the PN.S ob- 
servations, the velocity is a function of the separation between 



the position of the PN in the two frames. In this way, the slit- 
less technique avoids the two-stage selection and measurement 
process. Because there are no slits or fibres, the light loss is re- 
duced. The number of measurable PN velocities is not limited 
by the number of available fibres or the restrictive geometry of 

the slits. 

NGC 1399 is a cD galaxy in the Fornax cluster dSchomberj 

lT986h . It is one of the closest cD galaxies and therefore promises 
one of the best kinematic data sets for the outer halo. Located 
near the center of their clusters, cD galaxies are believed to form 
through accretion and mergers with other cluster members. As 
such, these galaxies are valuable clues to the formation history 
of ga laxies and the evolution of clusters dHausman & Ostrikeri 
Il978h . 

This paper describes the latest kinematic data for the outer 
parts of the galaxy. It is part of a larger project intended to dy- 
namically model the galaxy using PNe and X-ray constraints. 
NGC 1399 is we ll studied with kin ematic measurements from 
long-slit spectra dSaglia et al ] l2000h . Previous to this work, the 
radia l extent of the PN kine matics was limited to the velocities 
from lArnaboldi etail d 19941) at 24kpc. The current sample mea- 
sures more than 4 times as many PN velocities out to more than 
twice the radius. The resulting kinematics provide a significant 
improvement in understanding the velocity field of the halo. 

We describe our observations in §2 before moving on to the 
data reduction, the image registration and the detection of plan- 
etary nebula in §3. We discuss the measurement of the PN ve- 
locities in §4, and §5 addresses the technique for mapping the 
astrometry on spectral images. The analysis of the sample, in- 
cluding discussions of objects associated with NGC 1404 and 
a low-velocity subcomponent, is covered in §6. We explore the 
kinematics of the halo of NGC 1399 in §7 and give a summary 
of this wo rk in §8. Th r ougho ut the paper, we assume the dis- 
tance from lTonrv et all (1200 lh of 19.95 Mpc, which means that 
1" corresponds to 96.7 pc. 



2. Observations 

2. 1 . Instrumental set-up 

Observations of five fields around NGC 1399 were carried 
out with FORS1 (lAppenzeller et al 1 119981) on the 8-m VLT at 
Paranal Observatory during the nights of 25-27 Nov 2000. The 
five fields represent a subset of the 3x3 grid proposed, which was 
centered on the galaxy as shown in Fig.Q] Information from the 
observing log is tabulated in Table [1] The FORS1 field-of-view 
is 6'. 8 x 6'. 8 or 39.5kpc on each side at 19.95 Mpc dTonrv et al.l 
1200 lh . The spatial scale is 0". 25 /pixel. 

We used the 600B grating with a measured mean dispersion 
of 1.19A pix -1 . The interference filter, FILT_503_5+86, is cen- 
tered near the wavelength of the redshifted [OII I1 emission at 
the sy stemic velocity of NGC 1399 (1425 km s'. lGraham et al.l 
1998). The filter, with a FWHM of 60A, limits the length of 
the spectra of continuum sources which would otherwise cover 
much of the field. It allows detection of PNe with velocities 
within +1800 km s _1 of the systemic velocity. For some fields, 
bright stars were masked out. 

With this instrumental configuration, the direction of dis- 
persion is along the horizontal axis (x-axis) depicted in Fig. [2] 
Because of the slitless technique, both axes have spatial infor- 
mation. In the focal plane, the presence of the grism causes an 
anamorphic distortion resulting in a contraction in the direction 
of dispersion (see discussion in 32.31 ). 
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Fig. 2. Part of the central field (field 5). The top frame is the WO exposure and the bottom frame is the counter-dispersed E180 
exposure after rotating back in the computer. The diffuse light at the top center of each frame is the southern tip of the bright 
central region of NGC 1399. In these images, the stars appear as streaks in the x-direction (the direction of dispersion). Examples of 
planetary nebula candidates are circled. Notice that there is a small difference in the PNs' x-positions in the WO and El 80 fields due 
to their velocity. The arrows point to emission-line objects that are associated with some continuum- these are background galaxies. 



Table 1. Observing log* 



Field 


P.A. Label 


Exp. Time 


Seeing 






(s) 


(FWHM) 


NGC 1399-2 


WO 


3 x 1200 


1.42" 


NGC 1399-2 


E180 


3 x 1200 


1.34" 


NGC 1399-4 


WO 


3 x 1200 


1.23" 


NGC 1399-4 


E180 


3 x 1200 


1.10" 


NGC 1399-5 


WO 


3 x 1200 


1.11" 


NGC 1399-5 


E180 


3 x 1200 


1.09" 


NGC 1399-7 


WO 


3 x 1200 


1.04" 


NGC 1399-7 


E180 


3 x 1200 


1.22" 


NGC 1399-8 


WO 


3 x 1200 


1.09" 


NGC 1399-8 


E180 


3 x 1200 


1.10" 



Each field is observed at P.A.=0° and 180° 
respectively. 



labeled W0 and El 80 



width is related to the seeing in y-direction and the linewidth 
and the seeing in the x- (wavelength) direction. The location of 
the PN on the spectral image is a function of its position on the 
sky, its velocity, and the instrumental distortion. See Fig. [2] for 
an example image. In the second exposure (El 80), the instru- 
mental set-up remains the same except that the spectrograph is 
rotated to P.A. = 180°. In this configuration, the sky appears ro- 
tated 180° on the plane of the CCD — in order to register the two 
images we must counter-rotate the El 80 frame during the soft- 
ware reduction. See Fig. [3] for an illustration of the technique. 
The velocity is related to the difference in position between the 
PN in W0 and the counter-rotated El 80 frame (E180 raf ), after 
the images have been registered using the stars as shown in the 
final frame of Fig. [3] 



2.3. Calibrations 



2.2. Observing mode 

In this technique, the positions and velocities are obtained using 
two narrow-band, dispersed exposures of the same field. In the 
first, an image is acquired at position angle (P.A.) 0° (W0). The 
result is that stars become streaks in the direction of dispersion 
with a length determined by the dispersion and the width of the 
filter. Emission-line sources, such as PNe, appear as dots whose 



The position of each PN in the dispersed image is a combina- 
tion of its position on the sky, its line-of-sight velocity and the 
distortion from the instrument; we must quantify and separate 
these contributions in order to measure the velocity. The grism 
introduces an anamorphic distortion into the field that needs to 
be corrected before we measure velocities. Also, the dispersion 
changes as a function of position and the filter bandpass is de- 
pendent on the location in the field. 
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Fig. 3. Schematic describing the counter-dispersed imaging technique. In the field shown, there are four stars (shown as rainbow 
continuum sources) and 3 planetary nebulae (shown as small emission-line dots). The first frame shows the WO field with dispersion 
to the right. The second frame, E180 ro/ , shows the 4 stars in the same locations but now dispersed in the opposite direction. This 
image was taken at P.A.= 180° and rotated back during the reduction. Notice that the PN change position between these frames 
because of their velocity. We combine the dispersed and counter-dispersed images in the final frame by aligning the stars. The 
separation (Ax) of each pair of PN detections depends on the PN's velocity. 



We used three sets of daytime calibration frames based on 
a 10 x 19 grid of slitlets imaged through the narrowband [OIII] 
filter. In the first set, the slitlets are imaged in undispersed white 
light. This gives us the positions of the slitlets in the grid. In 
the second set, the white light is dispersed. This set determines 
the filter bandpass as a function of position on the CCD. The 
final set is illuminated with a dispersed HeAr lamp. There are 
three lines covered by our bandpass: /15015.67, /15047.74 and 
/15085.82. This frame is used to evaluate the anamorphic dis- 
tortion and to measure the local dispersion. 

The daytime calibration exposures for the following steps 
were taken over three days where each slitlet position was ex- 
posed once in each configuration (image of the slitlet, dispersed 
white-light spectrum, and dispersed HeAr spectrum). The slitlets 
are formed by 19 pairs of moveable fingers, each calibrated at 
10 positions across the field. All of the positions were co-added 
to create a master white-light image, a master white-light spec- 
trum, and a master HeAr lamp spectrum. The three master cal- 
ibration frames are shown in Fig. [4] The image processing and 
analysis necessary for these calibrations was done in IRAF and 
Mathematica. The y-position is taken to be the vertical center of 
each slitlet. To measure the position in the dispersion direction, 
40 rows were summed along the slitlet and a Gaussian was fit — 
the x-position is taken to be the mean of the Gaussian. See 94.4.21 
for a discussion of how well it is possible to measure the central 
wavelength of the filter bandpass. 

2.3.1 . Mapping the anamorphic distortion 

The first calibration step is to use the x,y-positions from the im- 
age of the slitlets in the first calibration frame and the positions 
of the brightest line (^5015.67) in the third calibration frame to 
map the anamorphic distortion introduced by the grism. With no 
distortion, the HeAr line would show a constant offset from the 
image of the slitlet on the CCD. We tabulate the difference in 
x and y and use it to make a map of the anamorphic distortion 
using the IRAF task geomap. A graphical representation of the 
matrix is shown in Fig. [5] All 190 slitlets shown were used in a 
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Fig. 5. Anamorphic distortion. Each black point represents the 
location of the midpoint of a slitlet in the undispersed image 
frame. The red points are the corresponding slitlets as measured 
using the /15015.67 line in the dispersed HeAr lamp. The entire 
red grid has been shifted by the mean offset in order to display 
higher order distortions rather than the rigid shift caused by the 
disperser. Notice that the grism has introduced a sort of pincush- 
ion effect with the largest effects being at the right and left edges 
towards the middle. 



5th-order polynomial fit. Using geotran, we applied the solu- 
tion from geomap to all of our images to remove the distortion. 

Effectively, this transformation converts the data from the 
spectral into t he image plane. We de fine the spectral and im- 
age planes like I Arnaboldi et al] d2007l) : the spectral plane is the 
dispersed image as observed, and the image plane has been cor- 
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Fig. 4. Calibration frames. Each vertical line is comprised of 19 nearly-adjacent slitlets. Using this 10x19 grid of slitlets, we imaged 
the grid (left), added a grism to the light path to disperse the white light (center) and illuminated the 190 slitlets with a HeAr lamp 
and dispersed the light (right). In combination, these three frames allow us to map the anamorphic distortion, measure the bandpass 
shift as a function of CCD position and calculate the local dispersion. 



rected for the distortion introduced by the grism thereby mak- 
ing an undispersed, distortion-corrected image for one wave- 
length (/15015.67). The PN identification and velocity measure- 
ment were done on the image-plane images because they are 
distortion-corrected. 

2.3.2. Calculating the local dispersion 

Next, we measure the distance between the 5015. 67A and 
5085. 82 A lines in the dispersed HeAr lamp image to derive 
the local dispersion. We assume that the dispersion is linear on 
this scale. The two-dimensional, 2nd-order function describing 
the dispersion at each of the 190 calibration slitlets is fit with 
surf it, and we interpolate to get a measurement of the local 
dispersion at any location. 

Because we measure the dispersion in the non-distortion- 
corrected (spectral) calibration frame, we multiply the local dis- 
persions by the linear term from the geomap solution in 32.3.11 
when we apply them to relate position to wavelength. This is 
possible because the higher-order terms are small compared to 
the linear term of the anamorphic -distortion fit. This converts 
our dispersions from the measured spectral-plane dispersions to 
image-plane dispersions. 

2.3.3. Measuring the bandpass shift 

Finally, we map the bandpass shift as a function of position. We 
do this by comparing the Gaussian center of the dispersed white 
light to the position of the /15015.67 line. We expect there to be 
an underlying smooth, 2-D function that describes the bandpass 
shift caused by the presence of an interference filter at a small 
inclination in a converging beam. We modeled this function with 
a 2-D cubic using Mathematica and now we can evaluate the 
bandpass at any location on the CCD. Fig.|6]shows the measured 
bandpass shift across the CCD in light blue. The smooth red 
surface represents the model for the underlying smooth function 
that determines the local bandpass. The rms of the residuals is 
less than 1 pixel. 

When we align the W0 and El 80 images using the stars, we 
must first correct for the bandpass shift. In the W0 and El 80 
images of the same star, the light path is different causing dif- 
ferent bandpass shifts in the two images. In general, the location 
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Fig. 6. Bandpass shift as a function of position on the CCD. 
The shift is measured by comparing the measured center of the 
bandpass to the location of the /15015. 67-line at 190 locations 
across the CCD. The bumpy, blue-toned surface represents the 
measured shift. The smooth red surface is a 2-D cubic fit to the 
blue data. We evaluate the fit when correcting for this shift in the 
stars. All three axes are measured in pixels. 



of a star in our images is measured with a 2-D Gaussian using 
the n2gaussfit task in IRAF. For each star used in the align- 
ment, we evaluate the two-dimensional cubic function at that lo- 
cation and apply the shift to the n2gauss£it coordinates. The 
bandpass shift is measured in the spectral plane, but we want to 
register image-plane images. Therefore, using geoxytran, we 
apply the anamorphic distortion correction from geomap to ob- 
tain the distortion-corrected bandpass-shifted location of each 
star. These new coordinates are used to obtain the rigid shift that 
will be applied to the distortion-corrected images before calcu- 
lating the velocity of the PNe. Because we measure the center 
of the filter relative to the /15015.67 line, it becomes our fiducial 
wavelength. The E 1 80 and W0 images will be aligned on the x, y 
positions of this wavelength in the stars, and the wavelength of 
the PNe emissions will be measured relative to that line. 
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3. Data reduction, registration and identification of 
emission-line sources 

3.1. Data Reduction 

The data reduction for the counter-dispersed technique was car- 
ried out in IRAF. Each field is imaged with three- 1200s ex- 
posures at each position angle. The frames are bias subtracted 
and flat fielded individually. The sky flats were taken at twilight 
using the narrow-band filter and the grism. The three frames 
at each P.A. were combined using imcombine with a standard 
sigma-clipping for cosmic -ray removal. 

The result is a total of 10 combined, bias-subtracted, fiat- 
fielded frames- 5 at P.A.=W0 and 5 at P.A.=E180. We subtract 
the background from each frame by smoothing with f median 
and subtracting the smoothed background from the spectral im- 
age. Next each of the frames is corrected for anamorphic distor- 
tion using IRAF task geotran as described in 32.3.11 - this trans- 
forms the frame from the spectral to the image plane. We rotate 
the El 80 image by 180° so that the stars are in comparable posi- 
tions in the W0 and E180 /o , frames; the calibration frames are no 
longer applicable to E180 roI because the direction of dispersion 
is opposite. 

3.2. Registering the images 

We must now register the W0 and E180 rof images in order to 
measure the positions and velocities of the PNe. The registration 
of the W0 and E180 ro/ frames is done with a rigid shift. Note 
that the shift is calculated from stars, which need to be adjusted 
for the bandpass shift as described in 32.3.31 More specifically, 
the frames are aligned on the virtual position of a set of ~ 10 
monochromatic stars at /15015.67. To achieve this, we measure 
the locations of the stars in the spectral plane, apply the band- 
pass shift, correct for anamorphic distortion for the shifted coor- 
dinates using geoxytran, rotate the distortion-corrected, shifted 
positions from El 80 only, and calculate the shift that will be ap- 
plied to the images. After the registration of the frames, we can 
create a difference image and identify emission-line objects. The 
resulting pair of registered images looks like Fig. [2] 

3.3. Identification of emission-line objects 

In these images, emission-line sources are easily distinguished 
from stars because they appear as points instead of streaks. For 
the identification we rely on the W0, E180, o/ and W0-E180, o , 
difference images of each field. Specifically, an object is added to 
the PN catalog if it is point-like in both the spatial and dispersion 
directions, and if it appears in both W0 and E180 ro/ at the same 
y-position and with similar intensity (it need not be exact since 
the seeing may vary between the two images). 

We identify three kinds of emission-line objects: 

- Those that are point-like in both wavelength and space- 
these are planetary nebula candidates; 

- Those emission-line sources that are associated with a 
continuum- these are background galaxy candidates includ- 
ing some Lyman-Q- galaxies at z=3.1. These objects can be 
either point-like or extended in the spatial direction; 

- Those with resolved or double emission in wavelength. They 
can be extended or point-like in the spatial direction. These 
are also background galaxies and include some [Oil] emit- 
ters. 





■ 






1 1 


1 1 



Fig. 7. Examples of the three classes of emission-line objects 
identified in our sample. Each frame is 1'.5 across and north 
is up. top: Four PN are circled. They are point-like in both di- 
rections, middle: The arrows point to two examples of back- 
ground galaxies. The top-right one is a continuum source that is 
extended in the spatial (vertical) direction. The bottom-left one 
shows an emission line superimposed on some continuum. This 
could be a Lyman-o- galaxy, bottom: The object in the box is 
slightly extended in the wavelength direction- we are seeing the 
double emission of an [Oil] emitter. 



Table 2. Breakdown of emission-line objects 



NGC 1399 


NGC 1404 


Low-vel. 


Unassigned 


Background 


PNe 


PNe 


PNe 


PNe 


Galaxies 


146 


23 


12 


6 


224 



Examples of these classes of emission-line sources are 
shown in Fig. [7] Further examples are available in Fig. [2] We 
identify a total of 411 emission-line objects in our fields. The 
breakdown of these detections into types of emission-line ob- 
jects is shown in Table [2] The definitions used to distinguish 
members of NGC 1399 from those objects belonging to NGC 
1404 and the low-velocity subcomponent are described in 34.61 
and 36.3l respectivelv. For the planetary nebulae, we can use the 
known rest wavelength of the emission, /15006.8, to acquire the 
velocities using measured separation between the W0 and El 80 
frames- this is described in 34.11 
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3.3.1. The efficiency of the technique 

We found a total of 187 PNe in the five fiel ds, including 146 
bound to NGC 1399. As a point of comparison. lTeodorescu et all 
d2005l) applied a related technique to NGC 1344, another Fornax 
galaxy, using t he VLT and found 1 97 PNe. In their technique 
(summarized in lMendez et alJ d200lh ). each field is imaged with 
broad- and narrow-band filters, and then a grism is added and 
the narrow-band image is dispersed. The velocity would be mea- 
sured from the separation between the PNe in the narrow-band 
image and the dispersed, narrow-band image. 

As elliptical galaxies in the same cluster, we expect the PN 
detection efficiencies in these two galaxies to scale with their 
color and luminosity. NGC 1399 is red compared to NGC 1344 
(B-V=0.95 versus B-V=0.87 ba sed on the extinction- c orrecte d 
total lumino sities published by Ide Vaucouleurs et"al ] (Il99im . 
According to Hui et al. I (fl993l the difference in color would in- 
dicate that the luminosity specific PN number density of NGC 
1399 is only 43% of that of NGC 1344. Also NGC 1399 is 
almost a full magnitude brighter which should have the ef- 
fect of nearly doubling the number of PN detected in NGC 
1399. Finally, Teodorescu et al. covered the area around NGC 
1344 uniformly, while the NGC 1399 observations only covered 
5/9ths of the proposed grid. We can estimate the number of de- 
tections we are missing by symmetrizing the detections from the 
fields we do have. This would increase our sample by about 55 
PNe. 

If we adjust our detection rate for the difference in color, 
magnitude and coverage, we would have detected 224 PNe, 
which is similar to the 197 NGC 1344-bound objects found by 
iTeodorescu et al.ld2005h . We conclude that the counter-dispersed 
spectroscopy technique described in this paper provides an ef- 
ficient strategy for collecting large PN velocity samples in the 
local volu me and is comparable to the singly-dispersed imaging 
technique dMendez et alj|2001b ITeodorescu et al.ll2005l) . 

4. Planetary nebula velocities 

4. 1 . Velocity calibration 

The data reduction and the registration were done using the po- 
sitions of the stars. As continuum sources, their spectra are ap- 
proximated by the nearly-symmetrical filter curve (see discus- 
sion in §14.4. 2\ and adjusted to compensate for the filter band- 
pass shift. Unlike stellar continuum sources, planetary nebulae 
are not subject to the filter bandpass shift and their positions on 
the CCD are dependent on both their positions in the sky and 
their velocities. The counter-dispersion causes each PN to ap- 
pear in different places in the WO and E180 ro/ relative to the 
stars (i.e. to the virtual position of monochromatic stars). This 
principle is illustrated in Fig. [3] The velocity is a function of the 
distance between the position of the PN in the WO and El 80 ex- 
posures and offset by the wavelength we selected for the stellar 
alignment (/15015.67). More specfically, 



where A is the measured wavelength of the planetary, Ao is a 
fiducial wavelength determined by the registration procedure (in 
this case, ,150 15. 67), dA/dx is the local dispersion and Ax is the 
separation, in pixels, between the PN locations in the WO and 
E180 fields. 

Using this calibration, we generate line-of-sight velocity 
measurements for the PNe in our sample. The mean measured 



velocity is different from the systemic velocity of the galaxy — 
this indicates a zero-point offset in the calibration which is 
discussed in 34.41 First, in 34.21 we will determine the errors. 
In 34.31 we seek to demonstrate that the velocities are inter- 
nally consistent and therefore represent an accurate 2-D relative- 
velocity field. 

4.2. Errors on the velocity measurements 

The individual velocity errors are a combination of the errors in 
reduction, registration and measurement. The total error is cal- 
culated using 17 duplicates- these are PN that are detected in 
the overlap region between two fields. For each set of two fields, 
the mean difference in the duplicates is estimated (see section 
14.51 ) and the individual velocity errors for our sample are deter- 
mined from the standard deviation of the differences between 
all the pairs from overlapping areas. We measure an rms in the 
difference of 52 km s _1 , corresponding to an individual velocity 
error of 37 km s _1 . Thi s value is comparable to the 40 km s er- 
ror on the velocities in iMendez et aiTd2001l) . but it is larger than 
the 20 km s" 1 errors in the Planetary Nebula Spectrograph data 
dDouglas et al.ll2002l) . 

4.3. Testing the internal consistency of the velocity 
calibration 

Due to the complexity of the calibrations for counter-dispersed 
imaging, we seek to confirm that our velocity calibration is cor- 
rect. The following three tests confirm that the velocity calibra- 
tion does not depend on the distance of the PN from the center of 
the galaxy, the x, y position of the PN on the CCD or the velocity 
of thePN. 

In the first test, we computed a running average of the veloc- 
ity as a function of radius for the PNe in the central field (field 5). 
Assuming point symmetry, we expect to see a flat mean-velocity 
curve since the running average relies on averaged annular bins. 
Since it is not physically possible for the annular bins to have 
a velocity gradient, a rising or falling curve would indicate that 
our velocity calibration depends incorrectly on radius. In Fig. [8] 
we observe a flat running average, which indicates that our ve- 
locity calibration is not dependent on the distance to the center 
of the galaxy. 

In the second test we looked for a trend in the velocity as 
a function of p osition on the CCD. We compared our velocities 
with those from lArnaboldi et alJd 19941) using 14 objects in com- 
mon between the two samples with a mean difference of 166 km 
s . We show the dependence in the difference as a function of x 
and y in Fig. [9] and we do not see a scaling with position. 

The final test is the measured difference between the mean 
velocity of NGC 1399 PNe and the mean of the NGC 1404 can- 
didates that we can see in field 7. If our velocity calibration were 
scaled with velocity, we would expect to see the difference in the 
means of the two galaxies scaled proportionately. Using a max- 
imum likelihood technique we fit the PNe in field 7 with two 
normal distributions, one for each galaxy. See Fig.fTUl 

We see a difference of 491+55 km s _1 between the peaks 
of NGC 1399 and NGC 1404. The fit shown in Fig. [TO] finds 
the peaks at 1260±47 km s _1 and 1751+28 km s . This sepa- 
ration is consistent w ith the separation of the values in the liter- 
ature of 522 km s _1 dGraham et al.lll998h . The spatial selection 
caused by the boundaries of field 7 may influence the mean ve- 
locity of the NGC 1404 objects detected. As discussed i n 36.21 
NGC 1404 displays rotation; according to iGraham et aUd 19981) 
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Fig. 8. Velocity distribution of our central field's PN candidates 
as a function of radius. The red line marks the running average 
computed with 20 PNe in an annular bin. The lack of a trend in 
the average indicates that the velocity calibration does not vary 
with radius. 
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Fig. 10. Histogram of velocities in field 7 is overplotted with 
the modeled velocity distribution. The model is a maximum- 
likelihood fit to the unbinned velocity distribution. We can com- 
pare the separation between the peaks to the difference in the 
published systemic velocities of NGC 1399 and NGC 1404. In 
this figure, we have applied the systematic shift of 166 km s _1 
discussed in 34.41 After the shift, the fitted galaxy velocities are 
1426 km s- 1 and 1917 km s _1 . 
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Fig. 9. 14 PNe in common between our sample and the 1994 
NTT sample as a function of position on the CCD. The triangles 
represent the difference in velocity as a function of x-position 
and the open circles show it as a function of y-position. The lack 
of a trend in either distribution indicates that the velocity calibra- 
tion is not a function of x or y on the CCD. The mean difference 
in velocity is 166 km s _1 with <x=88 km s _1 . These values are 
addressed in 34.41 



it would have the effect of lowering the velocities that we de- 
tect in NGC 1404 which would account for the slightly smaller 
measured separation. Our conclusion is that the velocity calibra- 
tion does not change as a function of distance to the center of 
the galaxy, position in the field or velocity. We conclude that we 
are accurately measuring relative velocities with this technique. 
We discuss the determination of the absolute velocity in the next 
section. 
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Fig. 11. Comparison of the current ve locity measurem e nts wi th 
those from the NTT measurements dArnaboldi et al.l d!994l) ). 
These 17 objects were detected in both samples. For the mid- 
dle range of NTT values, we see a consistent agreement with the 
more recent data. At the fast and slow extremes, the agreement 
deviates. In this figure, we have already applied an offset of 166 
km s , which is the mean offset in the 14 PNe at the center of 
the distribution (shown with the dashed line). 



4.4. Calculating absolute velocities 

As described in 34.31 the measured instrumental velocities are 
relative velocities. Absolute velocities are necessary to compare 
these measurements with other kinematic tracers. The mean ve- 
locity of the current sample is 1197 km s . Using a friendless 
outlier algorithm, we remove 12 low-velocity planetary nebulae 
that we believe are not virialized and therefore do not trace the 
mass of NGC 1399 (See 36.31 1. The mean instrumental velocity 
for the 146 PNe associated with NGC 1399 is 1250+17 km s" 1 . 

We compare the VLT counter-dispersed-technique velocities 
with the 1994 NTT m ulti-object-spectroscopy measurements 
dArnaboldi et al.|[l994l) . As shown in Fig. [TT] a comparison in- 
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dicates that the offset is constant for objects in the center of the 
NTT velocity distribution. Relying on the 14 out of 17 objects 
that were detected in both samples in this range, we calculate 
a systematic offset of 166 km s _1 . The dispersion in the differ- 
ence, 88 km s , is consistent with the combination in quadrature 
of the errors in the VLT and NTT samples: i.e. 37 2 + 72 2 88 2 . 
As the distribution in Fig. [9] illustrates, the offset between the 
previous work and this velocity calculation is not dependent on 
position on the CCD. 

The error in the offset is calculated from the dispersion, 
88/ Vl4 = 24km s _1 . When we apply the shift with this error 
to our measured values, we find a mean velocity of 1416+29km 
s 1 for our sample of PNe in NGC 1399. This is consistent with 
iGraham et al] dl998b where v jys =1425+9 km s 1 (NED) and 
iRichtleret ail d2004l) where v sys = 1441km s -1 for the GCs. Since 
the offset is constant, our measurements stand alone as relative 
velocities. For absolute velocities, we rely upon the adjustment 
of 166 km s _1 based on the NTT values. 



4.4.1 . Comparison to X-Shooter 

X-Shooter is a new, multi-wavelent h, medium-resol ution, 
single-target spectrograph for the VLT ( lD'Odoricol l2008). It is 
designed to cover a wavelength range of /13000-24800 using 
three arms: UVB (A3000-5595), VIS (^5595-10240) and NIR 
(/110240-24800). As part of the commissioning observations, we 
received observations of one of our PNe from the UVB and VIS 
arms. We were able to resolve the double [OIII] emission at 
/14958.9 and /15006.8. The long wavelength range and a mul- 
titude of skylines gives us a wavelength calibration with 1 A pre- 
cision and yielding a measured velocity for this PN of 1 1 17 ± 60 
km s . Having accounted for the systematic offset of 166 km s _1 
discussed in §14.41 this value is consistent with the measurement 
of 1059 + 37 km s _I from the counter-dispersed images. 



Fit to Filter 

Fit to Star with Filte 

Filter 

Star seen with Filte 




4950 



5000 5050 
Wavelength [A] 



5100 



Fig. 12. Comparison of the filter to a stellar template as seen 
through the filter. The filter response curve is shown as the dot- 
ted line and fit with a dashed Gaussian. The solid, gray line rep- 
resents a 5500K [Fe/H]=0 stellar template (IMunari et alJl2005h 
as seen through the filter. Notice the asymmetry caused by metal 
lines in the blue end of the spectrum. Comparing the dash-dotted 
Gaussian stellar fit to the dashed fit to the filter, we see that the 
measured position of the center will not be significantly affected 
by the absorption (< 0.3 A). 



4.4.2. Understanding the velocity offset 

In 34.3l we demonstrated that the velocity calibration is internally 
self-consistent. Now we d iscuss the systematic o ffset. Relative 
to the measurements from lArnaboldi et"ai1 (1 19941) . we see a dif- 
ference of 166 km s _1 . We explored several possibilities for the 
origin of this offset. 

The first possibility questioned the assumption that the short 
stellar spectra are symmetrical in wavelength. Depending on the 
temperature of the star, it seemed possible that metal lines were 
creating an asymmetrical feature on the blue side of the stellar 
spectrum as seen through the filter. Our assumption of symme- 
try is important because we are using the locations of stars in W0 
and E180, o/ to register the images. If asymmetry caused a sys- 
tematic shift of the center of the filter bandpass by say 2 pixels, 
it would introduce an error in the mean velocity of 131km s _1 . 

This issue was resolved by comparing the measured center 
of the filter to the measured center of a stellar template point- 
by-point multiplied with the filter bandpass. The centers were 
obtained by fitting Gaussians. The stellar templates are ATLAS9 
model atmospheres with l og(^)-4.5, no rotation and solar metal- 
licity dMunari et al.ll2005l) . We covered a range of 4000-8000K. 
A typical fit to the filter and the filter-template combination is 
shown in Fig. Q~2] We conclude that, over this range of temper- 
atures, the fitted center of the stellar spectrum compared to the 
Gaussian fit to the filter differs by less than 0.3 A or 20 km s 
and therefore the absorption lines are not the cause of the zero- 
point offset. 



The second possibility for the cause of the velocity offset is 
the influence of temperature on the filter bandpass. Interference 
filters in convergent beams generally experience a shif t of 1A 
to the blue for every 5°C drop in ambient temperature dJacobvl 
[1981 . Since the calibration frames are taken during the day, they 
could have introduced a systematic shift. However, according 
to the condition report, the dome was climate controlled to the 
night temperature during the day. We conclude that the temper- 
ature did not affect the velocity measurement. 

After completing the astrometry and confirming its accuracy 
to 0".3 in the dispersion direction (see §0, we recalculated the 
velocities using celestial coordinates. This approach cuts out the 
registration of the images since the astrometry is completed in- 
dependently for the two frames. Additionally, the stellar astrom- 
etry confirms that we understand the geometry of each plane 
(W0 and El 80) independently. For each PN, we obtained a, 6 
(J2000) in the W0 and El 80 frames. After calculating the sep- 
aration in arcseconds, we used the plate scale (0.25"/pix) from 
the astrometry to convert to a separation in pixels and subse- 
quently Angstroms using the local dispersion (~1.2 A/pix). The 
technique did not improve the agreement between our measure- 
ments and the 1994 NTT values. 

In conclusion, we have examined several sources for the ve- 
locity offset and none of them are viable. We are unable to ac- 
count for this difference between our measurements and the liter- 
ature. Because we are confident that our velocities are internally 
consistent, we approach them as accurate relative velocities. 
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Table 3. Duplicated detections from regions of overlap between 
fields 



Field 


Field 


# of Duplicates 


Velocity offset" 


2 


5 


3 


18kms~' 


4 


5 


10 


-24 km s" 1 


7 


8 


3 


-55km s _I 


4 


7 


1 


22 km s _I 


8 


phantom PNe 


16 


47 km s" 1 



The velocity offset is the amount necessary to add to the velocities 
from the first column in order to make them consistent with the 
velocities from the field in the second column 



4.5. Mosaicking the individual pointings on a global velocity 
field 

Some of the fields have regions of overlap, and there are 17 PNe 
detected in these regions. We use these duplicated detections to 
try to understand if the systematic velocity offset is global or if it 
changes field to field. We compared velocity measurements from 
each overlap region to determine the offset necessary to make the 
two fields in question consistent. 

With the central field (field 5) defined as our benchmark, we 
adjusted two fields (field 2 and field 4) using PNe duplicated 
in field 5. We applied a high-sigma criterion for rejecting the 
outliers. Field 2, the topmost field, contained only 3 usable PN 
duplicates for this purpose, and we calculated an offset of 18 km 
s between the objects in field 2 and field 5. Field 4, just to the 
east of the central field, contains 10 duplicates, and we calculated 
an offset of -24 km s _1 between the objects in field 4 and field 5. 
We apply the calculated difference to each field creating a global 
velocity field from the three local ones. Since these adjustments 
are on the order of the error in the velocity measurement, we find 
these three fields consistent with a global systematic offset. The 
duplicates are summarized in Table [3] 

The two southernmost fields do not contain reliable overlap 
with the fields above (see Fig.Q]). We created a smoothed, sym- 
metrized velocity field using the now-consistent velocities from 
the three known fields (see 87. II regarding the construction of 
this velocity fields). By assuming point symmetry of the PNe in 
fields 2, 4, and 5 and smoothing, we create a velocity field that 
covers the southern pointing (field 8). The field was evaluated at 
the location of 16 observed field 8 PN candidates thereby cre- 
ating an equal number of phantom PNe. We then adjusted the 
mean velocity of the field 8 candidates to match the mean ve- 
locities of our phantom PNe by adding 47 km s -1 . Now, having 
made the field 8 velocities consistent with the rest of the field, 
we used 3 duplicate PNe to calibrate the velocities of the field 7 
candidates, adding 55 km s _1 to bring them onto the same plane 
as fields 2, 4, 5, and 8. There is one PN in common between field 
4 and field 7, but we favored the more robust adjustment to field 
8 instead. This procedure ensures that the velocity measurements 
in all of our fields are self-consistent. 

4.6. The separation of NGC 1399 and NGC 1404 

As shown in Fig. Q~| the central part of NGC 1404 is observed 
in the southeast field (field 7). The PNe in this field come from 
both NGC 1399 and NGC 1404. The published mean velocities 
of these galaxies differ by over 500 km s , but, due to the dis- 
persions, it isn't possible to separate the PNe unambiguously by 
their velocities. We have adopted a probability of membership 



method that assigns each PN to a galaxy only if it has 95% prob- 
ability of belonging. 

The first criterion of membership is the relative surface 
brightness contribution from the two galaxies at the location of 
the PN. The light profile of NGC 1399 is taken from Dirsch et al 
(2003). The light profile for NGC 1404 is an r 1/4 -law fit to pho- 
tometry data from the same paper. The location of the planetary 
was then projected onto the major axis of each galaxy and both 
light profiles were evaluated at the respective radii. The ratio of 
flux at the projected distance of the PN is the ratio of probability 
of membership based on light alone. 

The second criterion is the relative likelihood that an object 
of a given velocity belongs to the velocity distribution of NGC 
1399 versus NGC 1404. The total velocity distribution of field 
7 was fitted with two Gaussians through a maximum likelihood 
technique (Fig. ITOb . Then, reversing the technique, the relative 
likelihood of membership for each velocity was calculated. The 
ratio of likelihoods provides a ratio of probability of membership 
from velocity only. 

Our two criteria were multiplied to produce a ratio of proba- 
bility of membership value. A planetary is assigned membership 
if it is 19 times more likely to belong to one galaxy than the 
other (95% confidence interval). The 95% confidence criterion 
results in only 6 objects of ambiguous membership. If we lower 
our standard to 60% confidence of membership we would have 2 
ambiguous PNe — 2 more PNe would belong to NGC 1404 and 
2 more to NGC 1399. We retain the more stringent criteria. In 
order to examine the specific cases of the ambiguous PNe, see 
the figures in j36.ll The Field 7 candidates contribute 9 members 
of NGC 1399, 23 members of NGC 1404, and 6 of ambiguous 
membership. This analysis results in a total of 146 NGC 1399 
planetary nebulae in our sample. 



5. Astrometry 

The astrometry was performed using the imcoord s package in 
IRAF . We visually matched stars from USNO-B dMonet et al.l 
120031) to those in our images (about 10 per field). Because they 
are spectral images, for each star we fitted a 2-dimensional 
Gaussian and shifted the location to account for the local band- 
pass shift before building up a matrix of star positions. It is not 
necessary to correct for the anamorphic distortion because that 
transformation is folded into the mapping from the spectral plane 
to the celestial plane. Our images are described with a 2nd-order 
solution relative to the USNO-B stars. We expect our absolute 
astrometry to be no b etter than the errors on the USNO-B coor- 
dinates which is 0".2 dMonet et al.ll2003l) . 

To check the validity of our astrometry we selected several 
stars that were not used in the astrometric fit and compared their 
calculated VLT (i.e. this sample) positions to the USNO-B co- 
ordinates. The residuals are <0".4. We also completed an as- 
trometric plate solution for a narrow-band NTT image covering 
roughly the same field. When we compare the NTT astrometry 
to the USNO coordinates, the residuals are 0".3 — comparable 
to those for the VLT data. This indicates that the spectral im- 
agery of the VLT data does not cause an unpredictable behavior 
in the geometry that would affect our velocity calibration. 

For each planetary nebula, we average the position in the W0 
and El 80, o , frames to get velocity-independent position coordi- 
nates. We applied the astrometric solution to these coordinates 
yielding an a, 6 (J2000) catalog for the emission-line objects. 
This catalog will be included in a subsequent paper. 
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Fig. 13. The locations of the PNe in space. PNe moving faster 
than the systemic velocity are marked with black crosses. Those 
moving slower than the systemic velocity of NGC 1399 are 
marked with black circles. The low-velocity objects are marked 
with solid red squares. Blue crosses indicate objects associated 
with NGC 1404 and green crosses mark the ambiguous PNe 
that are not assigned to either galaxy. North is up and the ma- 
jor axis lies at 110° (dotted line). The d otted ellipses repres ent 
the C=25mag/arcsec 2 isophotes based on lDirsch et al.l d2003l) . 

5. 1 . Comparison to X-Shooter 

As discussed in 34.4. II a PN from our sample was included in the 
commissioning data for the UVB and VIS arms of X-Shooter. 
The pointing coordinates for the PN were taken from the as- 
trometric solution from the counter-dispersed images. The slit 
placement required a blind offset because PNe would not be vis- 
ible in an acquisition exposure. As expected, the PN lands 1/3 
of the way along the 10"-long slit. It appears very bright, so 
we assume that the object is well-centered on the slit. The X- 
Shooter observations indicate that our absolute astrometric solu- 
tion is near the limit imposed by the USNO-B coordinates. 

6. The PN velocity sample 

6. 1 . Phase-space diagrams 

In this section we view the PN 2-D velocity field in various 
projections in phase space in order to investigate the kinematic 
properties of the galaxy. In particular, we present our PN sam- 
ple in four subcomponents: those associated with NGC 1399, 
those associated with NGC 1404, those ambiguous ones that 
have not been assigned to a galaxy and a low-velocity subcom- 
ponent which we believe is not bound to either galaxy. The NGC 
1404 and ambiguous candidates are defined in 34.6l and the low- 
velocity subcomponent is defined in 36.31 The spatial locations 
of the planetary nebulae are shown in Fig. [13] Notice that the po- 
sitions are primarily limited by the locations of the five pointings 
shown in Fig.Q] 

Comparing the line-of-sight velocity to the projected radius 
in Fig. [14] we see the main NGC 1399 population in black with 
the measured systemic velocity of this subcomponent marked 
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Fig. 14. Phase space distribution. The NGC 1399 planetary neb- 
ulae velocities are shown as a function of radius in black. The 
12 low-velocity objects are marked in red. The open blue cir- 
cles show the NGC 1404 PNe and the green squares are those 
of ambiguous membership (see 34.61 ). The systemic velocities of 
the two galaxies are marked in dashed (NGC 1399) and dotted 
(NGC 1404) lines. The position of NGC 1404 is marked with a 
solid blue line at r=588". Notice that a PN can be assigned to 
NGC 1399 based on its PA. despite being at a radius and veloc- 
ity consistent with NGC 1404. 

with the black dashed line. Notice the paucity of PNe at 300" — 
this is caused by the spatial coverage of the survey. There is a 
low-velocity subcomponent moving at about 700 km s slower 
than the galaxy. Referring back to Fig. [13] these objects are dis- 
tributed across the north-central regions of our observations. As 
confirmed in Fig. [8] the average velocity of the main population 
is constant at all radii. 

6.2. PN kinematics of NGC 1404 

The PNe assigned to NGC 1404 (marked in blue in Fig. [13]) dom- 
inate the sample in the southwest. In Fig. [14] these objects are 
found at large radii, around the syst emic velocity of NGC 1404 
(1947 km s ': lGraham et al.l (119981) ). As illustrated in Fig. [10] 
the fitted mean vel ocity for the NGC 14 04 PNe in this sample is 
1917 ±37 km s~'. lGraham et al] (Il998l) report rotation in NGC 
1404 that could cause the PNe that we are detecting to be up 
to 100 km s _1 slower than the systemic velocity of NGC 1404. 
This may contribute to the ~30 km s -1 difference between the fit 
shown in Fig.[10]and the mean velocity value in the literature. 

6.3. Low-velocity subpopulation 

We selected the 12 low-velocity objects from our sample us- 
ing a high-sigma cl ipping algorithm similar to the technique in 
iMerrett et al. (120061) . It is an iterative outlier test that rejects ob- 
jects that lie > 2.75cr from the mean of the total sample. The 
clipping limit of 2.75cr is chosen because 99.3% of a normal 
distribution lies within 2.75cr. In a sample of 146 objects, we 
expect to find less than one object outside this range. This sub- 
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population of 12 objects has a mean velocity of 718 km s . 
Removing these low-velocity outliers symmetrizes the sample 
around NGC 1399's v jyj because there is no corresponding high- 
velocity counterpart. 

In the high-sigma clipping algorithm, we assign PNe mem- 
bership in either NGC 1399 or the low-velocity subcomponent. 
It is possible or even likely that one or two of the objects as- 
signed to the subcomponent are actually members of NGC 1 399. 
However, in most of the analysis in this paper, it is important to 
define the sample as a discrete set of objects, and the results will 
not change significantly based on the misclassification of one or 
two objects. This is not the case for the calculation of the ve- 
locity dispersion, which is sensitive to outliers in the wings of 
the distribution. The velocity dispersion calculation employs a 
robust fitting technique which is detailed in 37.2.11 

The low-velocity objects are remarkable in that they are not 
uniformly distributed over the entire galaxy, but rather concen- 
trated in the north-central parts. Considering the asymmetric 
placement of the fields and the relatively small numbers of PNe 
in some of the outer fields, we considered the possibility of a 
selection from the filter. Because it is an interference filter at an 
angle in a converging beam, the bandpass function is not cen- 
tered on the CCD, as shown in Fig. [6] We examined the 50% 
transmission edges of the filter as a function of distance from 
the intersection of the plane of the CCD with the axis of the 
filter, as seen in Fig. [15] We conclude that the filter is not caus- 
ing the asymmetry in the velocity distribution. The most recent 
globular cluster (GC) work also shows a handful of low veloc- 
ity objec ts in the sample remo ved with a tracer-mass-estimator 
method (ISchuberth et al.ll2009l) . which supports our claim that 
we are seeing a real substructure. 

A low-velocity subcomponent is difficult to detect with in- 
tegrated light measurements. Based on the phase-space distribu- 
tion of the low-velocity PNe, we can infer that the subcomponent 
would manifest itself as an asymmetric feature in the line profile. 
This retrograde tail would be visible in the h 3 value from inte- 
grated light measurements. ISaglia et al.l d2000l) present h3 values 
of about -0. 1 at large radii along the minor-axis slit. These mea- 
surements agree with the asymmetry in the PN distribution tend- 
ing towards low velocities, but they do not provide an estimate 
for the amount of light or the mean velocity of the subcompo- 
nent. 

The origin of the low-velocity subcomponent is unknown. 
The mean velocity is si milar to the systemic ve locity of NGC 
1396 (808+22 km s -1 , iDrinkwater etail (1200 ll) ). an elliptical 
galaxy to the west of NGC 1399, and their spatial distribution 
is consistent with their association with the outer halo of this 
galaxy. Unfortunately, our observations do not extend to the 
west, so we lack more conclusive evidence. The low-velocity 
globular clusters are dominated by blue metal-poor members. 
This is consistent with the possibility of the low-velocity ma- 
terial originating from a stripped galaxy. We do not detect a 
velocity gradient in the low-velocity PNe, so it seems unlikely 
that these objects are being currently stripped as part of a in- 
teraction between the two g alaxies. The ph otometry of NGC 
1399 from the OBEY survey dTal et all2009h indicates a smooth 
light distribution down to the surface brightness limit of V=27.7 
mag arcsec~ 2 , so the stellar component accompanying the low- 
velocity PNe must be extremely diffuse across our fields. This 
corroborates the indication from the velocity gradient that we 
are not seeing a clearly defined stream of PNe from an infalling 
galaxy. 

In any case, the asymmetry in the velocity distribution in- 
dicates that the low-velocity PNe are not currently phase-mixed 
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Fig. 15. The detection limits imposed by the filter. This figure 
shows all of the PNe detected in our five fields. The black circles 
represent the ones associated with NGC 1399, the low-velocity 
subcomponent and the unassigned PNe. NGC 1404 PNe are 
marked with a blue open circ le. The dotted line m arks the sys- 
temic velocity of NGC 1399 dGrahamet al.lll998h . The red and 
blue squares trace the red and blue 50% transmission limits of 
the filter respectively. At small distances, the mean wavelength is 
near the mean wavelength of the galaxy. As the bandpass shifts 
at the edges of the CCD, the sample gets closer to the edge of 
the filter. Since there are no detections at or near the edge of the 
bandpass, we conclude the filter is not blocking the high- velocity 
counterpart to the low-velocity subpopulation. Moreover, we are 
detecting NGC 1404 PNe at these high velocities — it would be 
possible to see NGC 1399 ones at the same wavelength if they 
existed. 



with NGC 1399. Like in M87 (iDohertv et al.ll2009l) . we see a 
superposition of kinematic components in the halo of this cen- 
tral cluster galaxy. This reflects the complexity of galaxy super- 
positions and interactions in cluster cores. The subcomponents 
could not have been disentangled without discrete velocity trac- 
ers; having PN velocities is therefore important for deriving a 
complete picture of the kinematics in the outer halos of bright- 
est cluster galaxies. The detection of subcomponents in galaxy 
halos is not restricted to cD galaxies. In cases like M31 or the 
MW, these structures can be traced with red giant branch (RGB) 
stars. In nearby clusters, we are unable to detect individual RGB 
stars, but it is possible to measure the velocities of individual 
PNe. Globular cluster velocities can also be used outside the lo- 
cal group, but we emphasize that PNe hold an important role 
because they are the only individual stars visible at several Mpc 
and they trace the stellar population directly. 

7. Kinematics of the halo of NGC 1399 

7.1. Rotation 

Now we consider the kinematic properties of the main NGC 
1399 sample. The smoothed, symmetrized velocity map shown 
in Fig.[l6]illustrates the rotation in the observed fields. It is made 
by reflecting the entire sample (defined as the PNe which are not 
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Fig. 16. NGC 1399 halo velocities. This smoothed, symmetrized 
velocity map shows the location of each PN with an x. The pho- 
tometric major axis is vertical in this configuration. The map 
indicates rotation of <K 50 km s . The kinematic major axis is 
offset from the photometric major axis- the kinematic axis twists 
at large radii. In this figure, the 12 low-velocity PN were not in- 
cluded since we assume that they are not part of the dynam ical 
system. The dotted ellipse represents 2R e dSaglia et alJl2000h . 



Table 4. Rotation of the Halo of NGC 1399. The photometric 
major axis is at PA=1 10°. 
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classified as members of NGC 1404, the low- velocity subcom- 
ponent or the ambiguous objects) in x, y and velocity and then 
smoothin g using an adaptive G aussian kernel smoothing as de- 
scribed in lCoccatoe t al. (2009|). The velocity field was fit with a 
cosine function in order to derive the direction and amplitude of 
rotation as a function of radius. The measured values for these 
quantities and their errors (computed by means of Monte Carlo 
simulations) are obta ined using the same process described in 
ICoccato e t al. ( 2009) and are listed in Table [4] We observe ro- 
tation of ~ 20km s along an axis offset from the photometric 
major axis (PA. = 110°). Additionally, the kinematic major axis 
displays a strong twist in the o uter regions similar to those de- 
scribed in lCoccato et all (120091) . 

This analysis is consistent with the low-amplitude of rota- 
tion (< 30 km s along the major axis) observed by Saglia 
et al. (2000) inside 1 '. In contrast to the rotation result from 
lArnaboldi et al.l (1 19941) . we see only a low amplitude of rotation 
in the outer parts — rotation of < 20 km s is a fraction of their 
value. A comparison to the globular cluster system's rotation can 
be found in 37.3.21 



7.2. The velocity dispersion 

Using the new PN velocities, we trace the dispersion out to 
large radii. PNe are particularly good for this task because they 
extend well beyond the integrated-stellar-light measurements, 
and their number de nsity follows the light where they overlap 
(ICoccato et al.ll2009h . 

7.2.1. Robust Fitting 

Due to possible contamination from the low-velocity subcompo- 
nent, the dispersion is calculated using a robust fitting technique. 
We begin by measuring the standard deviation of the data. To 
avoid using the tails of the distribution, we estimate the disper- 
sion using only objects within 2 standard deviations (2cr) and 
then scale the measured dispersion using a factor generated in 
a Monte Carlo simulation of the same process. A 2<x limit is 
selected because we are confident that the contamination is not 
significant within this range. A step-by-step description of the 
process follows in the next two paragraphs. 

Based on the number of PNe in each annular bin, 53 objects 
are drawn at random from a Gaussian distribution with the mea- 
sured mean and dispersion of the data. Objects at greater than 2<x 
are rejected and the dispersion of the remaining objects is mea- 
sured. This process is repeated 10,000 times creating a Monte 
Carlo distribution of dispersion values. The mean of the 10,000 
dispersion values is used to calculate a scaling, fi, between the 
dispersion as measured within 2<x (0-2) and the known input dis- 
persion {cr lrue ). We define 



@~true 

The dispersion of the 10,000 measurements gives the error in our 
dispersion values, cr(cr). 

Returning to the data, a 2<x limit is applied to each radial 
bin and the dispersion is calculated for the PNe within this limit. 
This dispersion is then scaled by the factor defined in Eqn. [2] 
Since we expect the initial estimate for 2cr to be influenced by 
the outliers, we repeat this process until the dispersion estimate 
stabilizes. The error, cr(cr), is /2X the dispersion of the Monte 
Carlo distribution. 

Fig.fTTIshows the line-of-sight velocity dispersion as a func- 
tion of projected radius calculated with the robust-fitting tech- 
nique. The calculated values are included in Table [5] The radial 
limits of the bins are set such that all the bins have about the 
same number of objects (53 or 54). We discuss the dispersion, <x, 
rather than v, m ! , but they are approximately equivalent since the 
rotation velocity is low. The velocity error has been subtracted 
in quadrature. 

7.2.2. Comparison to previous PN measurements 

In co ntrast to the previous measurements (lArnaboldi et al.l 
1 19941) , the velocity dispersion from the new observations is much 
lower and much flatter. There could be several contributing fac- 
tors to this change. Compared to the 1994 values, this sample 
is four times larger, has more precise measurements and does 
not have a problem with the extreme velocities as described in 
Fig. QT| Moreover, the robust fitting treatment has a significant 
effect on the dispersion. This is shown as the difference between 
the red squares and the red line in Fig. [17] 

In summary, this reassessment of the velocity dispersion im- 
plies a lower value than the previously published measurements. 
In the latest NGC 1399 data, there is no evidence of an upturn in 
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Fig. 17. Dispersion profile of NGC 1399. The new data are shown as red squares in the outer regions. The instrumental dispersion 
has been subtracted in quadrature. For comparison, we have shown the calculation of th e dispersion s witho ut using the robust 
fitting (red dashed line). The black circles represent integrated light measurements from ISaglia et"al ] (120001) (solid =maior axis, 
hoilow =rninor axis). The minor axis is projected for comparison to the major axis. The PN measurements from lArnaboldi et"aD 
(fl99l are in green, rising steeply. The outermost integrated light point, a blue triangle, is from P. Saha (private communication), 
and it shows good agreement with a profile extrapolated from Saglia et al's (2000 ) major axis dispersion measurement and the new 
PN kinematics. The black lines are various models from lKronawitter et alJ d2000h (see discussion in text). 



Table 5. New velocity dispersion measurements in the outer halo 
of NGC 1399. The radius of each bin is the mean of the radii for 
all the objects in that bin. 
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1 115 
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Revised GCs based on data from (Schuberth et al. 2009) 
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224 


32 


2 242 


215 


31 


3 364 


230 


34 


4 508 


214 


31 


5 792 


176 


26 


Integrated Light from P. Saha (private comm.) 


163 


184 


+27/-8 



the dispersion which would have indicated a growing influence 
from the cluster potential. 



7.3. Other tracers 
7.3.1. Integrated Light 

The current PN dis persion values agre e well with the major-axis 
values described in Sag lia et al.l(l200 0). One side of the parallel- 
to-minor-axis values does not agree well with the PN measure- 
ments, which suggests that there may be preferential contamina- 
tion from non-phase-mixed PNe in this measurement. See Fig. 



The new PN measurements agree well with the outermost 
dispersion measurement from integrated light — less than 200 
km s _1 at 163"(P. Saha, private communication). This value, 
marked with a blue triangle in Fig. [TTJ is derived from one all- 
night, carefully-sky-subtracted, long-slit observation taken with 
the RGO spectrograph on the Anglo-Australian Telescope. This 
measurement is included in Table [5] Approximately equal time 
was spent on the target and the sky. The current PN sample is in 
good agreement with the largest-radius integrated light measure- 
ment confirming the consistency of these tracers. 

7.3.2. Globular Clusters 

NGC 1 3 99 has a rich glob ular cluster system with about 7000 
objects dDirsch et al.l 12003b . Their abundance makes them po- 
tentially valuable tracers, but only the red subsample follows 
th e radial distribu t ion of the light and the PNe (see Fig. 15 
in lSchuberfh et al.l (120091) ). We consider the velocity dispersion 
from the red clusters only — the blue clusters have a higher dis- 
persion which goes hand-in-hand with their flatter spatial distri- 
bution. Although the entire GC system is moving in the same 
potential as the PNe and stars, the different spatial distribution 
of the clusters requires different kinematics if they are to be in 
equilibrium. 

The GC dispersions included in ISchuberth et al.1 (120091) are 
higher than the PNe and Saha's integrated light value. These 
values are shown with magenta crosses in Fig. [18]The differ- 
ence seems to be highly influenced by contamination from NGC 
1404. The decontamination in that paper is based on a cut in 
radius — any objects within 3' of NGC 1404 were not included. 

We apply the two-fold decontamination al gorithm (as de- 
scribed in ^4.61 > to the red GC sample (RI in ISchuberth et alJ 



McNeil et al: Counter-dispersed spectroscopy of NGC 1399 



15 



500 



FT 



'co 400 h 



5 300 - 



g 200 



100 



i r 





J L 







200 



400 

R(") 



600 



800 



Fig. 18. Globular cluster dispersion profile of NGC 1 399. The black dots represent the same integrate d light measurem e nts de scribed 
in Fig. [T7] The magenta crosses connected with a dotted line mark the GC dispersion values from lSchuberfh et"aj] (l2009t) . After 
applying the decontamination algorithm from ^4.6l to Schuberth et al.'s RI sample and also removing two low-velocity GCs based on 
the tracer-mass estimator method in that paper, we measure the dispersions shown in magenta squares. An instrumental dispersion 
of 70km s" 1 has been s ubtracted in quadrature. This dispersion profile is in excellent agreement with one of the models from 



iKronawitter et al.l (120001) . as di scussed in §|7.4 | Linger ing contamination from NGC 1404 seems to have strongly influenced the 
velocity dispersion profile from lSchuberth et all (|2009i) around 250". 



( 120091) ) which leads to a sa mple of NGC 1399 GC s that has 
10 fewer objects than RIII in ISchuberth et al.l (120091). al t hough 
some of the objects that were rejected bv lSchuberth et alJ (120091) 
were kept in our sample. We also reject two low-velocity GCs 
from the RI sa mple based on the trac er-mass estimator (TME) 
algorithm from lSchuberfh et al.l(l2009l) . The GCs are rebinned to 
have an equal number of objects in each bin, but the dispersion 
values fall with in one error bar i f we us e the bin limits described 
in Table A.l of I Schuberth et al.l (120091) . After the decontamina- 
tion, the TME outlier rejection and the subsequent rebinning, the 
measured dispersions drop to ~ 200km s _1 and are in excellent 
agreement with the PNe. See Fig. [19] The revised GC velocity 
dispersion values based on the two-fold decontamination are in- 
cluded in Table |5] 



The globular c l uster rotation analysis is summarized in §7 
of ISchuberth et al.l (120091) . In agreement with the PNe, they find 
marginal rotation for the red sample. The red globular cluster 
rotation is 61+35 km s at an azimuthal angle of 154+33 de- 
grees, but a rotation of this amplitude could be measured from 
randomized data with probabilty 27% — this is interpreted as 
insignificant rotation. They find rotation of 110+53 km s _1 at 
®o =130°+24° in the blue sample, but we do not expect the blue 
GCs to agree with the rotation of the PNe and the integrated light 
based on t he spatial distribution a nd kinematics shown in Figs. 
14 & 15 of ISchuberth et al.1 (120091) . While there is disagreement 
about the position angle, the planetary nebula and globular clus- 
ter kinematics agree that the amplitude of any rotation in NGC 
1399 would be small. 



7.3.3. Fornax cluster galaxies 

Just as PNe and globular clusters can be used to trace the mass 
distribution of a galaxy, cluster members allow us to estimate the 
cluster's mass. NGC 1399's position at the core of the Fornax 
cluster leaves open two possibilities for the dispersion profile 
at extremely large radii; the halo of a cluster's central galaxy 
could be dynamically separate from the cluster, or it could be 
part of it. For example, in the case of the Virgo cluster core, M87 
displays a truncation at 160 kpc, wh ere the velocity disp ersion 
of the galaxy drops below 80 km s _1 (iDohertv et al.1120091) while 
the potential determined from the X-rays is already dominated 
by the cluster mass distribution, implying that the central galaxy 
is dynamically separate from the cluster. If a central galaxy is not 
truncated, then it may be difficult to define the boundary between 
the galaxy and the cluster since the velocity dispersion of the 
galaxy could rise with radius to become similar to that of the 
cluster. 



Photometry of NGC 1399 from ISchombertl (1 19861) indicates 
that there is no truncation, at least in the stellar mass, up to 
a radius of about 300kpc. The component of the Fornax clus- 
ter centered on NGC 1399 has a dispersion of 370 km s 
dDrinkwater et al.ll200ll) . As discussed by Drinkwater et al., the 
dwarf galaxies in Fornax constitute a distinct population with a 
much higher dispersion than the giants (429± 41 km s _1 com- 
pared to 308+30 km s _1 ). Fig. [T9l shows the veloci ty disper- 
sion f or likely Fornax Cluster members taken from iFergusonl 
(fl989l These values agree very well with Drinkwater et al. 
The most distant dispersion measurements of NGC 1399 avail- 
able c ome from the GCs- t he profile seems to be falling at 
1000"( ISchuberthetal.l l2 009l). At -2000" the dispersion pro- 
file is not well understood. Bergond et al.l (120071) see an indica- 
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Fig. 19. Dispersion profile representing the best values for each 
type of tracer. As in Fig. [17] the new PNe data are shown 
as red squares and the black ci r cles r epresent integrated-light 
measurements from lSaglia et al.l (120001) (solid=major axis, hol- 
low=minor axis). The last integrated-light measurement (blue 
triangle) agrees well with the new PNe dispersion s. For com- 
parison, the red GC sample from ISchuberth et al.l d2009h has 
been decontaminated with the same technique as the PNe and 
is displayed with magenta aster isks. The black asterisks repre- 
sent Fornax Cluster dispersions (lFergusonlll989l). The solid line 
shows the best model for the PN data from iKronawitter et alj 
(120001) . 



tion of the onset of t he clu ster potential dBergond et al.ll2007l) . 
but ISchuberth et al.l (120081) interpret the same data as being 
part of the NGC 1399 system with a falling dispersion profile. 
Current observations do not bridge the space between the last 
GCs (~2000") and the first dispersion estimate from the clus- 
ter members at nearly 3000". At the radial limit of our sample 
(500"), the PN dispersion reaches 220 km s _1 indicating that 
these PNe are not in the regime dominated by the cluster po- 
tential, but still within the NGC 1399 halo. Mass models based 
on AS CA and Rosat X-ray data dlkebe et alJll996hlPaolillo et al. 
120021) indicate that the transition to the cluster potential might be 
at around 500". However, these mass mo dels may not be consis - 
tent with the globular cluster kinematics (ISchuberth et alj|2009l) : 
this issue requires more work. 

7.4. Dynamical Models 

A full dynamical model based on the new PN and GC data 
is deferred to a subsequent paper. Here we only compare 
the da ta briefly to a range of models from IKronawitter et alj 
J2000I) . These models were constructed based on absorption 
line kinematics reaching 100", assuming spherical symmetry. 
The anisotropy of the orbit distribution was constrained by line 
profile shape measurements inside 80". Best-fitting distribution 
functions were constructed in a sequence of potentials including 
that of the stars with constant mass-to-light ratio and an addi- 
tional logarithmic halo potential. Best fit composite potentials 
and confidence regions were determined from the combined fit 




Fig. 20. top: C i rcular velocity curve for the model from 
IKronawitter et alj (HOOO) that best fits the inner absorption line 
kinematics, the PN velocity dispersions and the red GC velocity 
dispersions (solid line). This corresponds to the model shown 
with a solid line in Figs. [17] and [18] The dotted line shows the 
circular velocity curve of a constant mass-to-light ratio model, 
bottom: B-band mass-to-light ratio profile for the best-fit model. 

to all the absorption line data. The lines overlayed on Fig.[T7lrep- 
resent several of the models from lKronawitter et al.1 (|2000); note 
that that their projected velocity-dispersion profiles are available 
only to R = 420". The long-dashed model at the top repre- 
sents the best fit to the photometry a nd integrated-light kinemat- 
ics described in ISaglia et"ai] d2000l) . At the bottom, the dotted 
line shows a self-consistent model with constant mass-to-light 
ratio. The short-dashed line near the middle shows an interme- 
diate model that is consistent with the inner integrated light data 
as well as the outer-most PN point. The best match to all of the 
velocity dispersion points from the current PN sample and also 
to the globular cluster points fo r R < 400" is provided by a 
model at the low-mass end of the IKronawitter etaD ([2000) con- 
fidence range, shown with a solid line in Figs. [TT] & [18] This 
line represents the velocity-dispersion profile of a mildly radi- 
ally anisotropic (J3 0.25) spherical model with a central B- 
band mass-to-light ratio of M/Lb — 12, a dark matter halo with 
circular velocity v c 340 km s for R > 200" and a mildly 
increasing M/Lg-profile, reaching M/Lg - 27 at R ^ 800" (see 
Fig. 120b. Outside 100" the anisotropy of this model is not nec- 
essarily a good approximation for the halo of NGC 1399, how- 
ever, because of the shallow luminosity density profile the ve- 
locity dispersions are insensitive to anisotropy dGerhardlll993l : 
IChurazovet"ai1l2010l) . Nonetheless, this model is only indica- 
tive; more detailed dynamical modeling is needed to explore the 
full range of anisotropy and mass profiles consistent with the 
new data. 



8. Summary and Conclusion 

Using counter-dispersed, slitless spectroscopy from FORS1 on 
the VLT, we detected 411 emission-line objects around NGC 
1399, including 146 PNe associated with NGC 1399, 23 asso- 
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ciated with NGC 1404, 12 low-velocity PNe, 6 unassigned PNe 
and 224 background galaxies, and we measured the line-of-sight 
velocities of the PNe. This technique, applied for the first time 
on an 8-m class telescope, measures individual velocities with 
errors of less than 40 km s _1 at nearly 20 Mpc distance. It pro- 
vides an efficient way of sampling the velocity field in the outer 
halos of galaxies, particularly because it does not require a two- 
stage detection-and-measurement process. In this paper we de- 
scribed the reduction technique specifically developed for the as- 
trometry and velocity calibration of these observations. The new 
PN velocity measurements allow us to examine the kinematics 
of the outer halo of the cD galaxy NGC 1399. 

The 146 PNe in NGC 1399 show weak rotation. The kine- 
matic major axis is offset from the photometric major axis and 
displays a twist at large radii. The small am plitude of rotation 
confirmed the result from Saglia et al. | d2000h . The PN velocity- 
dispersion profile flattens at about 220 km s extending to 
~ 40kpc. The PNe dispersions show good agreement with the 
maj or-axis integrated-light values and the globular cluster sam- 
ple dSchuberth et al.1 l2009h after the NGC 1404 objects have 
been removed using a different method. 

A com parison to the sequence of dynamical models con- 
structed by iRronawitter et ail d2000h from inner absorption line 
kinematics indicates that the PNe in the outer regions of NGC 
1399 are consistent with a model with a B-band central mass-to- 
light ratio of 12 and a dark matter halo with a circular velocity of 
v c ^ 340 km s . This comparison highlights the importance of 
having kinematics at large radii to distinguish between the vari- 
ous models. This paper establishes the most up-to-date data set 
with integrated light, PNe and GCs, which will constrain new 
dynamical models to be constructed in a subsequent paper. 

The PN kinematics of the halo support the role of cD galax- 
ies in hierarchical structure formation within clusters. The multi- 
component sample is comprised of objects associated with NGC 
1399, NGC 1404 and a low-velocity sub-group at -800 km s -1 . 
At the cluster core, we expect to find a heterogenous popula- 
tion including stars left over from previous accretion, merger 
or tidal stripping events; these stars are still in the process of 
phase-mixing. The superposition of the various components on 
the NGC 1399 field underscores the complexity of velocity mea- 
surements in cluster cores and the necessity of using discrete 
tracers to detect these components. 

These results have opened valuable opportunities for further 
research on the formation and evolution of elliptical galaxies and 
on the dynamical structure of NGC 1399 and the Fornax clus- 
ter core in particular. The ambiguity of the source of the low- 
velocity data underlines the importance of understanding the 
kinematics of the western region of NGC 1399 and the implica- 
tions on the dynamics of the Fornax cluster core. Given that the 
outer regions of early-type galaxies are dark-matter dominated 
and the orbits in this range have long dynamical time scales, the 
PN kinematics will provide valuable constraints on dynamical 
studies of this galaxy. The new data indicates that at large radii 
we must be aware of contamination in the wings of the velocity 
distribution as we look towards modeling the galaxy. Using these 
new PN data and the existing GC and integrated light observa- 
tions, the dynamical modeling of NGC 1399 will be described 
in a future paper. Modeling NGC 1399 is a particularly excit- 
ing prospect because of existing X-ray data providing an inde- 
pendent assessment of the potential. The estimate of the orbital 
anisotropy will shed light on the formation of the halos of bright- 
est cluster galaxies. Finally, given the success of this technique, 
we will apply the developed reduction processes to other sys- 



tems to further explore the kinematics of the halos of early-type 
galaxies using PNe as discrete tracers. 
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